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LONG-TERM  GOALS 


The  goal  of  this  project  was  to  develop  an  automated  detection  and  tracking  algorithm  for  broadband 
targets  using  complex  hydrophone  data  from  a  passive  acoustic  array.  The  algorithm  is  an  integral  part 
of  a  larger  Coherent  Automated  Multi-Target  Tracker  (CAMTT)  system  under  development  by 
Metron,  Inc.  [1]  for  Detection,  Classification,  and  Localization  (DCL)  for  passive  Anti-Submarine 
Warfare  (ASW).  The  algorithm  integrates  the  Maximum  a  Posteriori  Penalty  Function  (MAP-PF) 
tracking  algorithm  [2]-[4]  with  the  Likelihood  Ratio  Detection  and  Tracking  (LRDT)  methodology 
[5],[6],  and  spectral  tracking  approach  of  [7],  The  detection  and  tracking  problem  is  treated  as  a  joint 
detection  and  estimation  problem  and  the  combined  system  automatically  (1)  detects  and  drops  targets, 
(2)  jointly  estimates  bearing  vs.  time  tracks  for  all  targets,  and  (3)  jointly  estimates  the  received 
spectrum  of  these  targets.  The  spectral  estimates  improve  the  detection  and  tracking  capability  and 
will  be  used  to  aid  the  classification  component  of  the  CAMTT  system. 

OBJECTIVES 

The  technical  objective  of  this  project  was  to  develop  an  automated  detection  and  tracking  algorithm 
for  broadband  targets  using  complex  hydrophone  data  from  a  passive  acoustic  array.  The  detection 
and  tracking  problem  was  treated  as  a  joint  detection  and  estimation  problem  and  the  combined  system 
automatically  (1)  detects  and  drops  targets,  (2)  jointly  estimates  bearing  vs.  time  tracks  for  all  targets, 
and  (3)  jointly  estimates  the  received  spectrum  of  these  targets.  The  goal  was  to  develop  an  algorithm 
that  outperforms  existing  techniques  in  terms  of  the  mean  square  error  (MSE)  in  tracking  detected 
targets,  increased  and  earlier  detection  of  quiet  targets,  and  reduction  in  false  detections. 

APPROACH 

The  detection  and  tracking  problem  was  treated  as  a  joint  detection  and  estimation  problem  operating 
directly  on  the  complex  hydrophone  data.  The  joint  detection  and  estimation  approach  is  theoretically 
difficult.  An  optimal  solution  is  intractable  to  derive  and  would  be  computationally  complex  to 
implement.  Traditional  multitarget  detection  and  tracking  techniques  therefore  partition  the  process 
into  two  isolated  processes:  broadband  target  detection  and  bearing  estimation  from  array  snapshot 
data,  followed  by  track  estimation  from  the  bearing  estimates.  This  partitioning  results  in  procedures 
which  are  suboptimal  and  require  data  association  to  match  bearing  estimates  to  targets.  More 
sophisticated  approaches  first  process  the  complex  hydrophone  data  to  provide  estimates  of  spectral 
content  across  frequency  and  beams  (frequency  azimuth  (FRAZ)  data)  and  use  the  beamformed  data  as 
input  to  the  tracker.  The  approach  is  still  partitioned  into  two  processes  and  must  still  incorporate  a 
data  association  component,  although  the  partition  occurs  closer  to  the  original  data.  In  order  to 
improve  performance,  a  new  approach  that  operates  directly  on  the  array  data  is  required. 

The  proposed  algorithm  integrates  the  MAP-PF  tracking  algorithm  [2]-[4]  with  the  LRDT 
methodology  [5],[6],  and  spectral  tracking  approach  of  [7],  The  algorithm  operates  directly  on  the 
complex  hydrophone  data. 

The  approach  of  MAP-PF  is  to  develop  a  likelihood  function  for  the  complex  hydrophone  data  as  a 
function  of  the  joint  target  state  of  all  detected  targets  and  then  to  find  the  joint  target  state  that 
maximizes  the  likelihood  function.  This  method  requires  knowledge  of  the  number  of  targets  and  an 
initial  estimate  of  the  bearing  and  spectrum  of  each  target.  A  key  feature  of  the  MAP-PF  approach  is 
the  use  of  the  penalty  function  method  of  nonlinear  programming  in  conjunction  with  an  alternating 
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maximization  procedure  to  obtain  a  tractable  solution.  Penalized  maximum  likelihood  estimates  of  the 
target  bearings  and  power  spectra  are  found,  which  then  act  as  a  synthetic  “measurements”  in  a 
Kalman  filter  update  of  the  state.  The  two-step  estimation  process  is  similar  to  traditional  methods, 
except  the  processes  are  coupled  via  the  penalty  function.  In  the  bearing/power  spectrum  estimation 
step,  the  current  target  state  is  used  to  guide  the  estimation  process  and  help  eliminate  ambiguous  or 
spurious  estimates,  and  to  eliminate  the  data  association  step  of  traditional  multitarget  tracking 
approaches.  In  the  track  estimation  step,  the  penalty  function  parameters  act  like  the  measurement 
error  variance  to  control  the  influence  of  the  bearing  and  power  spectrum  estimates  on  the  final  track 
estimate. 

The  LRDT  approach  is  used  to  determine  the  presence  or  absence  of  targets  and  provide  an  initial  state 
estimate  of  newly  detected  targets  to  the  MAP-PF  tracker.  The  LRDT  methodology  develops  a 
likelihood  ratio  between  the  density  of  the  complex  hydrophone  data  given  the  joint  target  state  of  all 
detected  targets  and  the  density  assuming  an  additional  hypothesized  target.  When  the  likelihood  ratio 
surpasses  a  threshold,  a  new  target  is  declared  present.  A  similar  procedure  is  used  to  determine  when 
an  existing  target  should  be  declared  absent. 

Combining  the  MAP-PF  and  LRDT  algorithms  is  straightforward  and  natural  since  they  both  compute 
the  same  joint  likelihood  functions.  This  approach  provides  a  complete  automated  algorithm  for 
detecting  and  tracking  targets  directly  from  complex  hydrophone  data,  avoiding  inefficiencies  inherent 
in  partitioned  methods. 


WORK  COMPLETED 

In  the  first  project  year,  we  had  the  following  technical  accomplishments: 

(1)  Developed  a  broadband  MAP-PF  tracker  for  multiple  targets  and  a  maneuvering  platform.  The 
target  state  is  the  bearing,  bearing  rate,  and  spectrum.  The  noise  background  is  assumed  known 
and  white. 

(2)  Developed  an  algorithm  for  computing  a  whitened  likelihood  function  required  for  the  LRDT 
detector.  The  algorithm  effectively  removes  known  targets  by  modeling  them  as  slightly  spatially 
spread  targets  rather  than  point  sources. 

In  the  second  project  year,  we  had  the  following  technical  accomplishments: 

(3)  Extended  the  MAP-PF  tracker  to  model  colored  noise  and  include  unknown  noise  parameters  in  the 
state  vector. 

(4)  Developed  a  version  of  the  algorithm  that  operates  on  beamformed  data.  This  provides  an 
intermediate  algorithm  that  separates  two  key  features  of  the  MAP-PF  algorithm 
underdevelopment:  the  MAP-PF  methodology  itself,  and  the  use  of  element  data. 

In  both  project  years,  we  also: 

(5)  Tested  the  algorithms  on  simulated  data  in  a  realistic  moving  target,  maneuvering  platform 
scenario. 

(6)  Assisted  Metron  in  integrating  the  algorithms  into  their  CAMTT  system  and  initial  testing  on  real 
data. 
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RESULTS 


An  effective  multi -target  tracker  for  broadband  targets  in  unknown  colored  noise  was  developed  and 
tested  on  simulated  data  in  a  realistic  scenario  involving  moving  surface  and  submerged  targets  and  a 
maneuvering  platform.  The  full  state  vector  consists  of  the  noise  parameters  and  a  target  state  vector 
for  each  target.  The  target  state  includes  the  bearing,  bearing  rate,  and  spectral  parameters.  Including 
the  noise  and  target  spectral  parameters  allows  the  tracker  to  discriminate  between  targets  that  are  close 
spatially  and  provides  a  significant  performance  improvement  for  tracking  “crossing”  targets. 

The  algorithms  were  tested  on  a  realistic  scenario  as  shown  in  Figure  1 .  There  are  three  targets:  two 
loud  surface  ships  and  one  quiet  submarine.  The  platform  is  a  surface  ship  towing  a  48-element  linear 
array  with  design  frequency  667  Hz.  The  targets  move  in  linear,  constant  velocity  tracks  in  3D-space, 
however,  tracking  is  done  in  bearing  space,  so  the  target  tracks  cross  and  appear  nonlinear.  The 
bearing,  bearing  rate,  and  total  target  power  as  seen  at  the  array  are  shown  in  Figure  2.  The  target 
spectra  vary  across  the  frequency  band  and  are  shown  in  Figure  3.  The  noise  consists  of  spatially 
isotropic  environmental  noise  that  decreases  in  power  across  frequency  plus  spatially  white  sensor 
noise  that  is  constant  across  frequency.  The  noise  spectra  are  shown  in  Figure  3.  The  array  collects 
broadband  acoustic  data  in  41  frequency  bins  from  200-600  Hz. 

During  the  bearing  estimation  step,  the  bearing  of  each  target  is  found  separately  by  maximizing  a 
penalized,  whitened  likelihood  function,  where  the  other  targets  are  considered  to  be  interference  and 
whitening  is  performed  with  respect  to  the  combined  noise  and  interference.  To  analyze  the  effect  of 
proper  noise  modeling,  we  pre-whitened  the  data  using  a  white  noise  model  fitted  to  the  data,  and  a 
combined  white/isotropic  model  fitted  to  the  data,  followed  by  conventional  beamforming  (CBF).  The 
result  is  shown  in  Figure  4.  The  thin  lines  are  the  CBF  output  in  each  frequency  bin  and  the  thick  lines 
are  the  CBF  outputs  combined  across  frequency.  It  can  be  seen  that  in  the  colored  noise  case,  the 
higher  frequency  likelihood  functions  contribute  more  to  the  combined  likelihood  function  because 
they  have  less  noise.  Also,  proper  colored  noise  modeling  allows  the  targets  to  stand  out  better  above 
the  noise  floor. 

The  whitened  likelihood  functions  for  each  target  used  in  the  MAP-PF  tracker ,  as  well  as  the  whitened 
likelihood  function  with  all  targets  removed  is  shown  in  Figure  5.  The  thin  lines  are  the  likelihood 
functions  in  each  frequency  bin.  The  thick  line  is  the  frequency-combined  likelihood  function,  and  the 
thick  dashed  line  is  the  peanlized  likelihood  function.  In  each  likelihood  function,  other  targets  are 
effectively  suppressed,  allowing  for  good  tracking  results.  The  final  whitened  likelihood  function, 
which  will  be  passed  to  the  LRDT  for  new  target  detection,  also  shows  effective  suppression  of  known 
targets.  The  target  tracks  in  bearing,  bearing  rate,  and  total  power  are  shown  in  Figure  6,  with 
excellent  results  comparable  to  the  results  obtained  last  year  with  white  noise  only. 

A  detailed  description  of  the  algorithm  will  be  published  at  the  IEEE  International  Conference  on 
Acoustics,  Speech,  and  Signal  Processing  in  March  2010  and  is  included  at  the  end  of  this  report. 
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3D  Scenario 


Figure  L  3D  scenario  with  three  moving  targets  and  a  maneuvering  platform . 


Figure  2 .  Bearing ,  bearing  ratef  and  total  power  of  targets  as  seen  at  the  array. 


Normalized  Spectra 


Frequency  (Hz) 


Figure  3.  Target  spectra  (normalized)  and  noise  spectra. 
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Bearing  (teg) 
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Figure  4.  CBF  output  with  (mismatched)  white  noise  whitening  and  colored  noise  whitening \ 


Figure  5.  Likelihood  functions  for  each  target  and  LRDT. 


Bearing  (deg)  Bearing  Rate  (deg/snapshot)  Total  Power  (dB) 


Figure  6.  Bearing,  bearing  rate,  and  total  power  tracks  overlaid  on  true  values. 
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Preliminary  results  for  MAP-PF  using  beamformed  data  show  that  it  does  not  perform  as  well  as  MAP- 
PF  using  element  data  during  target  crossings.  The  beamforming  operation  removes  phase  information 
that  is  critical  for  effective  whitening  of  other  targets  when  the  targets  are  close  together. 

IMPACT/APPLICATIONS 

The  work  is  directly  aimed  at  improving  the  performance  of  Navy  ASW  systems.  By  using  a  joint 
detection  and  estimation  approach  operating  directly  on  the  complex  hydrophone  data,  inefficiencies 
inherent  in  currently  used  partitioned  methods  can  be  avoided,  resulting  in  improved  tracking  of 
detected  targets,  earlier  detection  of  quiet  targets,  and  a  reduction  in  false  detections.  The  algorithm  is 
being  integrated  into  Metron's  CAMTT  system,  which  will  have  transition  capability  into  the 
Integrated  Undersea  Surveillance  System  (IUSS)  Integrated  Common  Processor  (ICP)  as  well  as  other 
systems  discussed  in  [1], 

TRANSITIONS 

The  MAP-PF  tracking  methodology  is  currently  being  implemented  for  tracking  frequency  in  the 
LandSafe  system,  a  laser  radar  system  for  determining  the  ground  speed  and  landing  distance  of  a 
helicopter,  which  is  under  development  by  Optical  Air  Data  Systems  (OADS),  under  an  SBIR  project 
sponsored  by  ONR. 

RELATED  PROJECTS 

GMU  is  not  currently  involved  in  any  related  projects. 

Dr.  Bell  was  acting  as  a  consultant  to  Optical  Air  Data  Systems  via  Signal  Processing  Consultants,  Inc. 
to  implement  the  MAP-PF  tracking  methodology  to  track  frequency  in  a  laser  radar  system  for 
estimating  helicopter  height  and  ground  speed. 
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ABSTRACT 

The  maximum  a  posteriori  penalty  function  (MAP-PF)  approach  is 
applied  to  tracking  the  bearing  and  bearing  rate  of  multiple  wide¬ 
band  sources  in  unknown  colored  noise.  The  track  estimation  prob¬ 
lem  is  formulated  directly  from  the  array  data  using  the  maximum 
a  posteriori  (MAP)  estimation  criterion.  The  penalty  function  (PF) 
method  of  nonlinear  programming  is  used  to  obtain  a  tractable  solu¬ 
tion.  A  sequential  update  procedure  is  developed  in  which  penalized 
maximum  likelihood  estimates  of  target  bearings,  spectra,  and  noise 
parameters  are  computed  and  then  used  as  synthetic  measurements 
in  a  set  of  Kalman  filter  trackers.  The  two  steps  are  coupled  via  the 
penalty  function.  During  parameter  estimation,  the  penalty  function 
prevents  erroneous  outlier  estimates  which  can  cause  the  tracker  to 
lose  track.  During  the  tracking  step,  it  determines  the  influence  of 
the  parameter  estimates  on  the  final  track  estimates  by  adaptively 
adjusting  the  measurement  error  variance.  Performance  is  demon¬ 
strated  on  a  typical  sonar  scenario. 

1.  INTRODUCTION 

Traditional  target  tracking  techniques  using  sensor  array  data  parti¬ 
tion  the  track  estimation  problem  into  two  isolated  processes:  bear¬ 
ing  estimation  based  on  data  collected  at  the  array,  followed  by  tar¬ 
get  tracking  using  the  bearing  estimates  as  measurements  [I].  The 
maximum  a  posteriori  penalty  function  (MAP-PF)  approach  [2]-[4] 
formulates  the  track  estimation  problem  directly  from  the  array  data 
using  the  maximum  a  posteriori  (MAP)  estimation  criterion  and  the 
penalty  function  (PF)  method  of  nonlinear  programming  to  obtain 
a  tractable  solution.  The  result  is  a  two-step  estimation  process  is 
similar  to  traditional  methods,  except  the  processes  are  coupled  via 
the  penalty  function  and  the  data  association  step  of  traditional  ap¬ 
proaches  is  eliminated  In  the  bearing  estimation  process,  the  penalty 
function  uses  the  current  target  states  to  guide  the  estimator.  In  the 
track  estimation  process,  the  penalty  function  determines  the  influ¬ 
ence  of  the  bearing  estimates  on  the  final  track  estimates  by  adap¬ 
tively  adjusting  the  measurement  error  variance. 

The  MAP-PF  tracking  techniques  in  [2]-[4]  were  developed  un¬ 
der  the  assumption  that  the  noise  was  spatially  white  with  a  known 
frequency  spectrum.  In  this  paper,  we  extend  the  technique  to  han¬ 
dle  the  case  of  unknown,  spatially  colored  noise.  In  order  to  focus  on 
the  noise  estimation,  we  assume  a  single  array  and  track  the  bearing 
and  bearing  rate  of  each  target  Extension  to  tracking  target  position 
in  three-dimensional  space  using  multiple  sensor  arrays  is  straight¬ 
forward  As  in  [4],  we  consider  wideband  sources  and  include  the 
target  spectrum  in  the  target  state. 


This  work  was  supported  by  the  Office  of  Naval  Research  under  Grant 
N000 14-08- 1-0225. 


The  statistical  model  is  presented  in  Section  11,  the  MAP-PF 
technique  for  unknown  colored  noise  is  developed  in  Section  111, 
simulation  results  are  presented  in  Section  IV,  and  conclusions  are 
given  in  Section  V. 

2.  STATISTICAL  MODEL  AND  ASSUMPTIONS 

The  model  consists  of  M  moving  targets  radiating  broadband  signals 
that  are  received  by  a  linear  array  of  N  sensors.  The  number  of 
targets  M  is  assumed  known  and  M  <  N.  Let  0  denote  the  target 
bearing  and  0  the  bearing  rate. 

At  the  array,  the  received  broadband  signal  is  sampled  and  trans¬ 
formed  into  the  frequency  domain.  We  assume  that  there  are  L  fre¬ 
quency  bins  of  interest.  Let  ~fkim>l  =  1, . . . ,  L  denote  the  received 
power  spectrum  of  the  mth  target  at  time  k.  The  (L+2)-dimensional 
state  vector  of  the  mth  target  at  time  k  is 

Skm—^Okm  Qkm  ~fk\m  ‘  “  ~fkLm  •  (1) 

We  assume  the  target  states  are  independent  of  each  other  and  the 
motion  of  the  targets  is  described  by  a  first  order  Gauss-Markov  pro¬ 
cess,  i.e., 

8km  =  FSfc_i,m  +  W  kmy  (2) 

where  w*m  is  a  zero-mean,  white  Gaussian  noise  process  with  co- 
variance  matrix  Q„,  and  F  is  the  block  diagonal  matrix 

F  =  block  diag  [Fi.It]  ,  (3) 

where  II  is  L  x  L  identity  matrix  and 


where  Af  is  the  time  interval  from  k  —  1  to  k.  Therefore,  the  proba¬ 
bility  density  function  (pdf)  of  Sfcm  given  8k  i  ?m  is 

P(S|fcm|st  im^A'fFSk  l,m,Q.)  (5) 

with 

p{ *Om  )  ~  A  (SomtOom)'  (6) 

The  N  x  1  observed  data  vector  in  the  lih  frequency  bin  during 
the  fcth  observation  snapshot  has  the  form 

M 

*kl  =  bktmVl(0km)  +  nkly  (7) 

m  =  l 

where  bkim  is  a  random  signal  sample  with  E[bktmbllfn]  =  7*1™. 
The  vector  vt(Okm)  is  the  N  x  1  array  response  vector  for  the  Ixh 
frequency  bin  to  a  target  whose  bearing  is  9km .  The  vector  n ki  is 


To  appear  in  tCASSP  2010. 


8 


an  N  x  1  vector  of  spatially  colored  noise  samples  with  covariance 
matrix  £[iu/n£{]  =  R(a *i),  The  noise  covariance  matrix  is  param¬ 
eterized  by  a  p-dimensional  parameter  vector  a whose  values  vary 
across  frequency  and  snapshots.  The  signals  and  noise  are  assumed 
to  be  sample  functions  of  stationary,  zero-mean,  Gaussian  random 
processes.  It  is  assumed  that  the  observations  are  collected  and  pro¬ 
cessed  so  that  they  are  independent  from  snapshot  to  snapshot  and 
frequency  bin  to  frequency  bin. 

To  handle  the  unknown  colored  noise,  we  define  a u  to  be  the 
p-dimensional  noise  state  vector  in  the  Ith  frequency  bin  at  time  k , 
i.e. 

a  ki  =  [am  *  *  •  «Hp]  •  (8) 

We  assume  the  noise  states  are  independent  across  frequency  and 
independent  of  the  target  states.  The  noise  state  equation  is  given  by 

(9) 

where  Vki  is  a  zero-mean,  white  Gaussian  noise  process  with  covari¬ 
ance  matrix  Qrt.  Therefore,  the  pdf  of  a  ki  given  a*  ltj  is 

p(&ki\&k-i,i)  ~ -V"(afc_u,  Qa)  (10) 

with 

p(aoi)  ~  N (ioi,  nSi) .  (11) 

At  time  k,  we  denote  the  collection  of  target  states  across  all 
targets  as  sir ,  and  the  collection  of  noise  states  across  all  frequency 
bins  as  a Given  s*  and  a k,  the  pdf  of  the  array  data  Zki  is  then 
jointly  complex  Gaussian, 

p(zki |s k ,a k)  ^  CM  (0, Kki (sfc , afc )) ,  (12) 


where 


Kfc/(sfc,afc)  =  7 klmVl(0km)yi(0km)H  +  R(a/fc*).  (13) 


Let  Z,  S,  and  A  denote  the  collections  across  snapshots  of  array 
data,  target  states,  and  noise  states,  respectively.  Assuming  there  are 
K  snapshots,  the  joint  pdf  of  Z,  S,  and  A  is  given  by 

p(Z,  S,  A)  -  p(Z|S,  A)p(S,  A) 

=  p(Z|S,  A)p(S)p(  A) 

=  (nn^i— ))• 

'  M 


( 

3,  MAP-PF  TRACKING  ALGORITHM 


(  K  \ 

l[p(Skm\Sk  l.m)lp(SOtn) 

n  p(«fci  |afc  p(aoi) 


(14) 


data.  First,  the  noise  parameters  are  assumed  known  and  penalized 
ML  estimates  of  the  source  bearings  and  power  spectra  are  obtained 
as  in  [4],  except  the  colored  noise  covariance  matrix  is  used  in  place 
of  the  white  noise  covariance  matrix  in  the  pre-whitening  step.  Next, 
noise  parameter  estimates  are  obtained  for  each  frequency  bin  as  in 
[7].  These  estimates  are  then  passed  to  the  KFs  as  measurements, 
and  the  target  and  noise  states  are  updated.  An  explicit  pseudo-code 
description  of  the  algorithm  is  provided  in  Table  1 . 

The  two-step  estimation  process  is  similar  to  traditional  meth¬ 
ods,  except  the  processes  are  coupled  via  the  penalty  function.  In  the 
parameter  estimation  step,  the  predicted  target  states  from  the  KF  are 
used  in  the  penalty  function  to  guide  the  estimation  process  as  well 
as  to  eliminate  the  data  association  step  [1]  of  traditional  multitaiget 
tracking  approaches.  In  the  track  estimation  step,  penalty  function 
parameters  act  like  the  measurement  error  variance  to  control  the 
influence  of  the  parameter  estimates  on  the  final  track  estimates 

The  penalty  function  parameters  are  chosen  to  be  equal  to  the 
Cramer-Rao  Bound  (CRB)  for  the  beanng/power  spectrum  estima¬ 
tion  problem  in  unknown  colored  noise.  The  CRB  will  be  smaller  for 
stronger  targets  than  weaker  targets,  and  will  be  larger  when  targets 
are  close  together.  This  provides  a  natural  tightening  or  relaxation  of 
the  penalty  function  for  the  different  target  and  noise  parameters  as 
the  scenario  evolves. 

Narrow  band  (single  frequency  bin)  Fisher  Information  (FI) 
terms  can  be  obtained  by  combining  the  analysis  in  [8,  pp.  961  -962] 
with  the  analysis  in  [7]  or  [9]  for  colored  noise.  The  bearing  FI 
components  are  then  summed  over  frequency  as  in  [10],  and  the 
wideband  CRB  is  found  by  inverting  the  resulting  joint  FI  matrix. 
The  CRB  C(s k ,  a* )  is  a  function  of  the  actual  target  and  noise  states 
s k  and  a*,  which  arc  the  quantities  being  estimated.  We  evaluate 
the  CRB  at  the  predicted  state  estimates  from  the  Kalman  filter  The 
penalty  parameters  are  then  obtained  from  the  appropriate  diagonal 
elements  of  the  estimated  wideband  CRB. 


4.  SIMULATION  RESULTS 

The  algorithm  was  tested  on  a  typical  sonar  scenario  in  which  a  mov¬ 
ing  platform  is  towing  an  N  —  48-element  linear  array.  The  array 
collects  broadband  acoustic  data  in  L  =  41  frequency  bins  from 
200-600  Hz.  There  are  two  high  signal-to-noise  ratio  (SNR)  targets 
and  one  low  SNR  target  The  array  platform  and  targets  move  in  lin¬ 
ear,  constant  velocity  tracks  in  3D-space,  and  the  platform  performs 
a  maneuver  and  changes  direction  during  the  observation  interval. 
A  bird’s  eye  view  of  the  scenario  in  xy-space  is  shown  in  Figure  1. 
In  bearing-space  where  the  tracking  is  performed,  the  target  tracks 
cross  and  appear  nonlinear,  and  move  across  array  endfire  during  the 
platform  maneuver  The  bearing,  bearing  rate,  and  total  target  power 
as  seen  at  the  array  for  each  of  the  targets  is  shown  in  Figure  2.  The 
target  spectra  vary  across  the  frequency  band  and  are  shown  in  Fig¬ 
ure  3.  The  noise  is  a  combination  of  isotropic  environmental  noise 
and  spatially  and  temporally  white  sensor  noise.  The  noise  covari¬ 
ance  matrix  in  the  l\h  frequency  bin  at  time  k  has  the  form 


We  wish  to  estimate  the  random  variables  S  (the  target  tracks)  and 
A  (the  noise  parameter  tracks)  from  the  observations  Z  We  find  the 
MAP  estimate  [3]  by  maximizing  p(Z,  S,  A)  with  respect  to  S  and 
A  .  Following  [2]-[4],  we  introduce  a  set  of  auxiliary  variables  n  and 
a  and  use  the  penalty  method  of  nonlinear  programming  to  make 
the  problem  tractable.  The  result  is  a  sequential  track  state  update 
procedure  based  on  the  Kalman  filter  (KF)  [6].  After  each  observa¬ 
tion,  target  and  noise  parameter  estimates  are  obtained  from  the  array 


R(aic/)  =  oVfcjR,,/  +  (13) 

where  R,  *  is  the  normalized,  known  isotropic  noise  covariance  ma¬ 
trix  in  the  Ith  frequency  bin.  The  noise  state  has  two  parameters 

a  ki  =  crl,,ki]  •  (16) 
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Table  1.  Sequential  MAP-PF  algorithm  pseudo  code. 


Initialize  so™  =  5om ,  P5|o,m  =  ;  ao/  =  ao/ ,  =  ft^ 

for  fc  =  1 , . . . ,  K, 

Predict  current  state 
for  m  =  1 , . . . ,  M 

—  1  ,m  =FSfc-l,m 

end  {m} 

for  /  =  1,. .  .,L 

&k\k  i,i  =afc*i,/ 

end  {/} 

Calculate  CRB  and  update  penalty  parameters 

C  k  =  C(sfc|fc-i,  ®fc|fc-i) 

e r«  =  [cs]m.m 

=  [(CZL,m](( 

sj/  =  [cs]lif 

for  m  =  1, . . . ,  M 
P  re -whiten  data 


for  /  =  1, . . . ,  L 

Tfe/m  =  |'>fc|fc-l,J,m'Vj(0fc|fc 

m'fkm  ( 

yf1  1  ,m')  +  R(afc|fc-u)} 

*«  =  rfc"lm2a« 

end  {/} 

Estimate  bearing  and  pow  er  spectrum 
L 

pkm  =  arg™X  ^{r?fc/m(^)  -lnWm(/i) 

I  /=1 

where  rjkim(p)  =  max 

and  vi(/i)  =  S“^2v,(/i) 
for/  =  1, . . .  ,L 

Pit/m  =  (£fcm)  —  1)  /  |v/(pjbm)|2 

end  {/} 

Pkm  —  [p Arm  Pfclm  ’  '  ’  PkLm\ 

end  {m} 

Estimate  noise  parameters 
for  /  =  1, . . . ,  L 


l,|*«v»(/i)|  /|v/(;j)|3 


=  "*£“*  I  (a,  A«,-  ‘  iMw/)  (from  [7]) 
end  {/} 

ACF  target  and  noise  state  updates 
H  =  [1  0  1  •  1]T 

form  =  1, . . . ,  A/ 

RSm  =  dtag([6ft.Ci -«£]) 

J>t|*-i,m“Fp;_1|*_limFT  +  Q. 

G;m  =  P'M*  1.mH7jHPJ|lt  lmHT  +  RJm}-* 
Sfcm  =  Sfc|fc-l,m  +Gfcm  {Ajkm  “  Hsfc|fc-l,m} 

Pa  _  r>a  pi  Tin' 

fc|fc,m  *  fc|fc— l,m  fc|fc  —  l,m 

end  {m} 

for  /  =  1, . . . ,  L 

Ra  _ 

kl  —  ** kl 

Rfc|fc-1,/  =  Rfc-llfc-l,!  +  Qa 

GS,  =  {PS|fc_lf,+R!,}  1 

a/t/  =  afc|fc-i,/  +  G£j  {ajt/  —  ajt|fc~i,/} 

Pa  _  na  pn  na 

k\k,l  —  rlt|fc-l,l“Wlrfc|fc-l,l 

end  {/} 
end  {fc} 


FIg.1.  Bird’s  eye  view  of  array  platform  and  target  trajectories. 


Fig.  2.  Bearing,  bearing  rate,  and  total  power  as  seen  at  the  array. 


Fig.  3.  Normalized  target  spectra. 
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Fig.  4.  Whitened  likelihood  functions  at  snapshot  fc=182. 


During  the  bearing  estimation  step,  the  bearing  of  each  target 
is  found  separately  by  maximizing  a  penalized,  whitened  likelihood 
function,  where  the  other  targets  are  considered  to  be  interference 
and  whitening  is  performed  with  respect  to  the  combined  noise  and 
interference.  The  whitened  likelihood  function  for  each  target  used 
in  the  MAP-PF  tracker  at  time  k  =  182  is  shown  in  Figure  4.  The 
dashed  line  is  the  penalized  likelihood  function,  where  values  away 
from  the  current  target  state  are  suppressed.  In  each  whitened  like¬ 
lihood  function,  other  targets  and  the  colored  noise  are  effectively 
suppressed,  and  good  bearing  estimates  are  obtained. 

The  bearing,  bearing  rate,  and  total  power  tracks  are  shown  in 
Figure  5,  with  excellent  results.  All  three  target  tracks  are  well  main¬ 
tained  through  target  crossings  and  when  the  targets  move  through 
array  endfire  during  the  maneuver.  Key  features  of  the  algorithm  that 
contribute  to  the  performance  are  using  the  CRB  as  the  penalty  func¬ 
tion  parameter,  and  including  the  target  spectra  and  noise  parameters 
along  with  the  DOA  parameters  in  the  state  vector  that  is  tracked. 


5,  SUMMARY 


Wc  have  developed  a  technique  for  tracking  the  bearing  and  bearing 
rate  of  multiple  broadband  targets  in  unknown  colored  noise  based 
on  the  MAP-PF  tracking  approach.  A  sequential  update  procedure 
was  developed  in  which  penalized  ML  estimates  of  target  bearings 
and  spectra,  and  noise  parameters  arc  found  and  then  used  as  syn¬ 
thetic  measurements  in  a  set  of  Kalman  filter  trackers  in  which  the 
target  spectra  and  noise  parameters  are  tracked  as  well  as  the  bear¬ 
ing  and  bearing  rate.  The  parameter  estimation  and  tracking  steps 
are  coupled  via  the  penalty  function,  which  is  critical  to  the  success 
of  the  technique.  During  parameter  estimation,  it  prevents  erroneous 
outlier  estimates  which  can  cause  the  tracker  to  lose  track.  During 
the  tracking  step,  it  determines  the  influence  of  the  parameter  esti¬ 
mates  on  the  final  track  estimates  by  adaptively  adjusting  the  mea¬ 
surement  error  variance. 


Fig.  5  Bearing,  bearing  rate  and  total  power  tracks  overlaid  on  true 
values. 
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